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Abstract 

We propose a system of nonlinear integral equations (NLIE) , which 
gives the free energy of the osp{l\2s) integrable spin chain at finite 
temperatures. In contrast with usual thermodynamic Bethe ansatz 
equations, our new NLIE contain only a finite number of unknown 
functions. On deriving NLIE, we use our osp{l\2s) version of the T- 
system and the quantum transfer matrix method. Based on our NLIE, 
we also calculate the high temperature expansion of the free energy 
and the specific heat. 
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1 Introduction 

In recent years, thermodynamics of solvable lattice models related to super- 
algebras have been studied by thermodynamic Bethe ansatz (TBA) method, 
and TBA equations are examined from various point of view (see, section 1 in 
for a comment on the present status of TBA equations for models related 
to superalgebras) . In our previous papers, we have derived TBA equations 
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for the osp{l\2) model 0, ^] and the osp{l\2s) model from the string hy- 
pothesis and the T-system ||^ of the quantum transfer matrix (QTM). These 
TBA equations contain an infinite number of unknown functions, and thus 
are not always easily treated. It is important to reduce the TBA equations to 
tractable integral equations which contain only a finite number of unknown 
functions. 

As for the XX Z spin chain, which is related to the algebra Ug{A^^^) of 
rank one, Takahashi proposed a nonlinear integral equation (NLIE) re- 
cently. This NLIE contains only one unknown function. Due to its simplicity, 
we can calculate the high temperature expansion of physical quantities 
to very high order from this NLIE. The purpose of this letter is to derive a 
system of NLIE with only a finite number (the number of rank s) of unknown 
functions from our osp{l\2s) version of the T-system f^. This letter is the 
first attempt to derive this type of NLIE for a vertex model associated with 
an algebra of arbitrary rank 

In section 2, we briefly mention the T-system for the osp{l\2s) model and 
the QTM method []T^ |T^, |TB|, 0. This section overlaps with section 2 



and section 4 in [y, but normalization of the fused QTM is different from 
In section 3, we derive the NLIE. The normalized QTM defined in section 2 
play the role of the unknown functions of these NLIE. These NLIE contain 
a parameter m, which corresponds to the fusion degree of the model. For 
m = 1, these NLIE form a closed set of equations, which contains only a 
finite number of unknown functions. On the other hand for m > 2, they 
couple with the ones for m = 1, and contain an infinite number of unknown 
functions. This type of equations (for m > 2) has never been considered be- 
fore even in the case of Takahashi's NLIE for the XX Z spin chain These 
NLIE relate to traditional TBA equations through the dependant variable 
transformation (|5.1|) . Using our new NLIE, we calculate the high temper- 
ature expansion of the free energy in section 4. Section 5 are devoted to 
discussions. 



^In the case of the deformation parameter q of an underhning quantum affine algebra 
is root of unity {\q\ — l,q ^ 1), the T-system becomes a finite set of difference equations. 
Thus the corresponding TBA equation becomes a finite set of integral equations. See, 
[|j for TBA analysis of integrable field theories related to osp{l\2s). We also note that 
different types of NLIE with finite numbers of unknown functions for different algebras of 
arbitrary rank were considered in 0, nOl in rather different contexts. 
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2 T-system and QTM method 



In this section, we consider [ffl] an integrable spin chain |jTT], O] associated 



with the fundamental representation of osp{l\2s), and introduce the QTM 



T3| , |T4|, [T^, |T6|, and the T-system for this modeL The i?-matrix[| 



19, 20, 21, 12] of this model is given as 



2v 

R{v)=I + vP-- E, (2.1) 

2v- g 

where g = 2s + l; v e_C; = i-lf'-''^P^''^5adSbc; = aab{a-^)^, a, 6, c, G 
{1,2, ...,s,0,s, ...,1}; l-<2^----<s^0^s-<---^2^1;ais 
(2s + 1) X (2s + 1) anti-diagonal matrix whose non-zero elements are aa;a = 1 
for a G {1, 2, . . . , s, 0} and aa,a = — 1 for a G {s, s — 1, . . . , 1}; a = a; p{a) = 
for a = 0; p(a) = 1 for a G {1, 2, . . . , s} U {s, . . . , 2, 1}. The Hamiltonian of 
the present model for the periodic boundary condition is given by 

H = J ^ iPk,k+i H — Ek,k+i \ , (2.2) 
fe=i ^ ^ ^ 

where J is a coupling constant: J > and J < correspond to the ferro- 
magnetic and antiferromagnetic regimes respectively; L is the number of the 
lattice sites; Pk,k+i and E^^k+i act nontrivially on the k th site and /c + 1 th 
site. 

The QTM is defined as 

JV 
2 

ti\v) = Trj Yl Ra^kjiu + iv)Ra^^_^,j{u - iv), (2.3) 

k=l 

where R'^{v) = Rl'^{v); Rjkiy) =^''Rkj{v) (tk is the transposition in the k- 
th space); N is the Trotter number and assumed to be even. By using the 
largest eigenvalue T]^^^(0) of the QTM ( |2.3| ) at f = 0, the free energy per site 
is expressed as 

/ = -T lim logT(')(0), (2.4) 

where we set u = (-^- temperature); the Boltzmann constant is set to 

1. One can obtain the eigenvalue formulae of the QTM ( ^.3] ) by replacing the 
vacuum part of the dressed vacuum form (DVF) for the row-to-row transfer 
matrix with that of the QTM. This DVF is imbedded into a DVF for a fusion 
hierarchy of the QTM. It reads as follows. 

m a 

Tln\v) = E n n v-'-{m-a-2j + 2k)), (2.5) 

{C(,fe}i = lfc = l 
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where the summation is taken over djk G {1, 2, . . . , s, 0, s, . . . , 2, 1} such that 
djk ^ c^j+i.fc and dj^ -< dj^+i- The functions {z{a;v)} are defined as 



z{a;v) = ipaiv) 



Qa.i{v + ^{a+l))Qa{v + ^ia-2)) 
Qa-iiv + ^ia-l))Qaiv + ^a) 

for a E {1,2, ... , s}, 



z(0; v) = My) ^; ■ \ , (2.6) 



Qa-ijv - |(a - 2s))Qaiv - |(a - 2g - 3)) 
Qa-i{v - |(a - 2s - 2))Qaiv - f (a - 2s - 1)) 
for a e {1,2, ... , s}. 



where Qa{v) = Ylk=i('" " '"fc"'')! ^ ^>o! Qo('y) := 1- The vacuum parts 
are given as follows 

/• 4>+(v)<t>- (v+i)'l>+ (v-^i) r ^ _ . 

U Mv--^i) ^ - ^' 

^a(t^) = <( CaMvH-iv) for 2^a^2, (2.7) 

,^_(.),^^(. 0^ (^+^0 for a=l, 

where (f)±{v) = {v ± iu)^; (a is a phase factor, {f^"^} is a solution of the 
Bethe ansatz equation (BAE) 



ad ) 



■>ad 



where fee {1,2,..., M^}; a G {1,2,..., s}; (5.+i(^^) := Qs(^^); = 2(5^^ - 
^a,d+i — ^a,d-\\ IS a phase factor. For a G {1, 2, . . . , s} and m G Z>i, we 
define a normalization function 



ATI") I 



2 VY'+V- I 2 

m a 



TT TT / m — a — 2j + 2k^ , m — a — 2j + 2k ^ 

xHllMv ^ ^)Mv ^ 0, (2.9) 

j=i k=i 

and set ft\v) = T^\v)/M^\v). For s = 1, f[^^ (f) coincides with eq.(4.8) 
in The poles of Tm\v) from the functions {Qb{v)} (dress part) are 
spurious under the BAE ( p.8|) . If one formally set the vacuum parts ( p.7|) 
and ( p^.9D to 1, then the remaining part (dress part) of the DVF Tm\v) is 
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same as the row-to-row case. Thus Tm\v) satisfies the following functional 
relation, which has essentially the same form as the osp{l\2s) T-system in 



t!^\v + -)Ti^\v --) = T;:i{v)Ti:Uv) + t^^-'\v)t^:-^'\v), 

for a e {1,2, ...,s- 1} and m G Z>i, (2.10) 



where 



T^^''\v) = 1 for aeZ>i, 
~rnw ^ (j)Jv - ^^i)(f)+(v + ^i) 

T^\^) = I ^^-7 ^ mGZ>i. 2.11 



3 Nonlinear integral equations 



In Takahashi's NLIE for the XXZ-model |^ was derived from the T- 



system of the QTM. Here we derive our new NLIE from our T-system ( p.lO| ). 



A numerical analysis for finite N, u, s indicates that a two-string solution 
(for every color) in the sector = Mi = Ms = ■ ■ ■ = M, of the BAE ( ^ 
provides the largest eigenvalue of the QTM ( p.3|) at f = [Q. From now on, 
we consider only this two-string solution. In this case, the phase factors are 
= Ca = 1 for any a. Moreover, we expect the following conjecture is valid 
for this two-string solution 

Conjecture 3.1. For small u (\u\ ^ 1), a G {l,2,...,s} and m G Z>i, 
every zero {zm^} ofTm\v) is located near the lines Imv = ±^2±Si^ _^ g+rn-a ^ 

fi^\v) has poles only at ±P^2\ (A; = 1,2): = - lu, j^t^ = 

^~^"^~"' i—iu. These poles are of order A^/2 at most. Moreover lim|i,|^oo ^m''('^) = 
Qm^ is a finite number 

nia) - ( {m + g)\m\ yr \ {k + a){k + g - a) \'' 

^"^ \{m + a)\{m + g - a)\ ) AJL \ k{k + g) J ' ^ ' 

This is a solution of the Q-system |l|: 

{Qt^f = Q'^:^_,Q'^:!^, + Qt'^Qt'-'^ for a G {1, 2, . . . , s - 1}, 

+ Qt-'v^:>, (3.2) 
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where m G Z>i; g[,"^ = Q^^) = 1. We are considering the case without an 



external field. If an external field exists, will be deformed. 
We may assume 

f 2 r Aa) -ria) >| 

where the coefficients b^m,j,ki'^^m.j,k ^ ^ given as follows 

^->=L + (3.4) 

Here the contour C^^^^ must surround counterclockwise manner and 

must not contain {p 7^ A;), "/^m,^) ~/^m,2j contour C^|^ must sur- 
round — counterclockwise manner and must not contain —^in!p {p 7^ k), 
/^mi5/^m2- Using the T-system ( p. 101) , we can modify ( |3.4|) as 



dv rxr!i(i;-|)TSi(^- 



where Tm^^^iv) := Tm\v). The first term and the second term in the first 



(resp. second) bracket {■ ■ ■ } in (p.5|) have common poles at Zm + i (resp. 
Zm^ —i)- However, these common poles are spurious since Tm\v) has no pole 
at these points. We also note that a subsidiary condition Ym\z^ ± |) = — 1 
for the excited state TEA equation (see, p. 2343 in and (|5.1|)) follows from 
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this observation. Substituting ( |3.5| ) into ( p.3| ), we obtain 




'^m-liy + f^ln!k 2 )^m+l ill + /^mi 2 ' 



for a G {1, 2, . . . , s} and m G Z>i, (3.6) 

where Tm^^\v) = Tm\v). Here the contour Cl^\ (resp. C^"*^) must surround 
counterclockwise manner and must not contain /3m,p— {p k), —j3^rni~ 

(resp. -~Ptl+~Pt\iP^k)rPt\+~^ti~^^^^^^^ For 
m > 2, the first term and the second term in the first bracket {■ ■ ■ } in ( |3.6| ) 
have a common singularity at which contributes to the contour integral. On 
the other hand for m = 1, this singularity at from the first term disappears 
since T^"'\y) = 1. Thus the contribution to the contour integral from the 
first term in the first bracket {■ ■ • } in (|3.6|) vanishes as long as the contour 
Cj"2 does not contain the singularities at zl""^ — (^fl + i (cf. Conjecture ^.11 ). 
This situation is parallel with the second bracket {■ ■ ■ } in ( p.6| ). Therefore 
for m = 1, (13.61) reduces to 



2 



k=l "^l.fc 



c(;)2vr. v-y-m T[''\y + ^ - 



dy ^ - [^J Tt'\y - pp + |)rr^^(i/ - /3S + 1) 

for a G {l,2,...,s}, (3.7) 
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where T^^^^\v) = Ti^\v). Here the contour C["'l (resp. C^il) must surround 
counterclockwise manner and must not contain z[^^ — ^[""l + i, — 
ip^k), - P^l - (resp. ~z[^^ + - ^, -/^g + P^l \p ^ k), 

P^Cl + PiiJ'Cl + P^i)- Now we shall take the Trotter limit iV ^ oo in (pj). 



k=l 



m,k 



- 

m,k 



k=l 



ri^\y + P':]k-^) 
dy 1 r 41 - /^i"), + |)rit (y - /^S + 



for a G {1, 2, . . . , s} and m G Z>i, (3i 
where %[^\v) := lim^^,. f 4"^); /^^ := lim;v-.oo /^i^ = /^^ ^ 

m- 
2~ 



4+'V) = 4V); 4"^(^) = 1; %i°\v) is a known function: 



4°^ (-) = Jim_ r(°) (.) = exp I - ^^."^i^^ I • (3.9) 



Here the contour (resp. C^"*;,) must surround counterclockwise manner 
with the condition I?/ 1 < \v—f3^j^\\ (resp. \y\ < and must not contain 

/3m,p — PU^i. {p ^ k), — — PU^i^, — /5^^2 ~ /^m,fc (resp. —I3m]p + {p ^ k), 
& + /^m,fc) /^m,2 + /^m,fc)- particular for m = 1, we have 



k=l 



rc[:l2m {v - y - f3[:l)T,'^\y + (3[^l - ^) 

+ dyTl^''\y-ptl + ^)rr'\y-p';l + i 



k=l 



for a G {l,2,...,s}, (3.10) 



where 4'*^^''(f) = Ti'^\v). Here the contour C["'l (resp. c["^) must surround 
counterclockwise manner with the condition \y\ < \v — f3[^l\ (resp. \y\ < 
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\v + I) and must not contain zi'^ - (3^1 + - f3'fl (p^k), -p['^l - 

= WmN^ooZ^i^- In ( CT ), O is written binomial coefficient: 



g(a) ^ ^2s+i^ ^ ^Q^g ^j^g f^^^ ^j^g^^ ^j^jg j^Ljg (|3jq ) contains only 

a /inite number of unknown functions {T^"'\v)}i<a<s- Moreover, we can 
express {T^\v)]i<a<s;m&>^ in terms of the solution of (|3.iq ) {T}''\v)] i<a<s 



by using a Jacobi-Trudi determinant formula 0. In this sense, we need not 
consider (|3.8|) for m G Z>2 in practical calculations. However, consideration 
on ( |3.8D for m G Z>2 manifests the relation between our new NLIE and the 
traditional TEA equations (see, section 5). The free energy per site is given 
as 

f = -J-T\ogT}'\^). (3.11) 



Solving ( |3.10 ) iteratively, we can obtain the free energy through ( |3.11| ). For 



s = 1 case, we find (|3.1CI|) converges numerically at least for \T / J\ > 0.37. 



4 High Temperature Expansion: osp{l\2) case 

It was pointed out |^ that Takahashi's NLIE for the XXX-model is very 
useful to calculate the high temperature expansion of the free energy. In this 
section, we shall calculate the high temperature expansion of the free energy 
( p.ll| ) for s = 1 by our new NLIE (|3.10| ) . We assume the following expansion 



for large T: 

rl'\v)=e^v{Y.<'n{v){iA. (4.1) 

\n=0 / 

Due to the shift of the argument y of Ti^\y) in (|3.10| ), we have to take into 
account the residues from the coefficients {a„(f)}, which contrasts with the 
XXX-model case [^. Thus the derivation is not so easy as the XXX-model 
case. But we still expect that it is easier than to use the traditional TBA 
equation which contains an infinite number of unknown functions (as for the 
high temperature expansion for the XXX-model from the traditional TBA 



equation, see for example, p. 124 in p9|). To calculate the coefficients {an{v)} 
efficiently, we need further assumption for n G Z>i: 

^ r b ■v'^i c 1 
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where bn,j,Cnj G C are independent of v. Substituting ( [4.1| ) into ( |3.10| ), we 
obtain the coefficients {a„(f)} up to order 12. For example, we have 

ao{v) = log 3, 

2 1 



ai{v) 



a2(v 



3iv^ + 1) 3(t;2 + f ; 



4(77 + 32i;2) 4(27 + 32^;2^ 



405(t;2+l)2 405(t;2 + |)2 



8(727 + 1040t;2 + 448t;^) 28 (837 + 720 1;^ + 128 1;^) 



10935(t;2 + 1)3 10935(^;2 + |)3 

-9097 + 781314 + 783969 + 175808 v^) 
7381125(f2 + 1)4 

144376263 + 138298104 1;^ + 37692144 + 2812928 1;^ 
59049000(t;2 + |)4 

Using {a„(0)}, we obtain 

/ , ^ 5 J 172 J2 20296 J3 52010 466964 
-log3 - — - - + + 



(4.3) 



T 27T 243 T2 59049 531441 1594323 

252697291 J6 2867981638 319036008559 



1937102445 T6 17433922005 1255242384360 
68608132529023 J9 69872025931694239 

^ 1921650566216724 T9 ^ 288247584932508600 T^o ^ ' '* 

157451559799196839 482908438587228085517^2 j 



648557066098144350 693437215072135939020 T12 ^ '^^^t' 

We also calculate the specific heat C = —Tg^. 

_ 344 J2 40592 208040 9339280 505394582 
~ 243 T2 ~ 19683 T3 ~ ITTMTT^ ^ 1594323 ~ 129140163 

40151742932 2233252059913 137216265058046 

+ 



13^ 



5811307335 156905298045 53379182394909 

69872025931694239 7^° 1731967157791165229 



+ 



3202750943694540 T^o 64855706609814435 T^i 
48290843858722808551 J12 ^ J ^3 

^ 5253312235394969235 T12 ^ '^'^T'' ^ ' 



We have plotted the high temperature expansion of the specific heat ( [4. 51) in 



figure 0. For comparison, we also plotted a pade approximation for (|3). For 



large T, ( [4.5|) agrees with our TBA analysis (see, Fig.l in 0]). This indicates 
the vahdity of our new NLIE ( PTTUQ . 
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Figure 1: Temperature dependence of the high temperature expansion of the 
specific heat for s = 1 and J = — 1: the broken hne denotes the plain series 
^31); the smooth hne denotes its pade approximation. 



5 Discussion 

In this letter, we have derived a system of NLIE with a finite number of 
unknown functions, which describes thermodynamics of the osp{l\2s) inte- 
grable spin chan. This type of NLIE for arbitrary rank is derived for the first 
time. We shall point out a relation between our new NLIE and usual TBA 
equations. The functions {Tm\v)} are related to the dependent variables 
{Ym\v)} of the F-system (or TBA equations) as 

^I'^H^) = 4%^^^^5?l^' = ~7^p4r;^^^ (5.1) 



for a G {1, 2, . . . , s — 1} and mGZ 



>i- 



This relation is also valid in the Trotter limit. 

In closing this letter, we shall enumerate future problems: 
• There are T-systems for other algebras [^, |2^, |2^, ^ |]. Using 



these T-systems, we can derive NLIE similar to the ones in this letter, which 



will be reported elsewhere BO 



• In this letter, we have considered the case without an external field. When 
an external field exists, a fugacity expansion of the free energy can be calcu- 
lated recursively through ( p.lO|) (as for the fugacity expansion of the XXX- 



model from the traditional TBA equation, see for example, p. 127 in 
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• One will be able to extend (|3.10| ) to trigonometric or elliptic case. In this 
case, one must take into account the periodicity in the summation in ( p.3|) . 
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